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ABSTRACT 

Aims. We study the nucleation of quark matter drops at the center of cold deleptonized neutron stars. These drops can be made 
up by unpaired quark matter or by color superconducting quark matter, depending on the details of the equations of state 
, for quark and hadronic matter. The nature of the nucleated phase is relevant in the determination of the critical mass M cr of 
hadronic stars above which it is possible a transition to a quark star (strange or hybrid). 

Methods. We investigate the dependence of M cr upon the parameters of the quark model (the Bag constant B, the pairing gap 
A, and the surface tension a of the quark-hadron interface) and for different parametrizations of the hadronic equation of state. 
We also calculate the energy released in the conversion of a pure hadronic star having the critical mass into a quark star. 
Results. In general, the dependence of M cr on B, A and a is mild if the parameters of the quark model correspond to hybrid 
stars, and strong if they correspond to strange stars. Also, the critical mass always decreases with A, and increases with B and 
a. The total released energy is in the range 3 x 10 52 erg - 4 x 10 53 erg. 

Conclusions. For a large part of the parameter space corresponding to hybrid stars, the critical mass is very close (but smaller 
than) the maximum mass of hadronic stars, and therefore compatible with a "mixed" population of compact stars (pure hadronic 
up to the critical mass and hybrid above the critical mass). For very large B the critical mass is never smaller than the maximum 
mass of hadronic stars, implying that quark stars cannot form through the here studied mechanism. The energy released in the 
conversion is sufficient to power a gamma ray burst. 
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1. Introduction 

The nucleation of quark matter in neutron stars has been 
studied by many authors, due to the potential connec- 
tion with explosive astrophysical events such as super- 
novae and gamma ray bursts. Some of the earlier stud- 
ies on quark matter nucleation (see e.g. Horvath et al. 
1992, Horvath 1994, Olesen & Madsen 1994, and refer- 
ences therein) dealt with thermal nucleation in hot and 
dense hadronic matter. In these studies, it was found that 
the prompt formation of a critical size drop of quark mat- 
ter via thermal activation is possible above a tempera- 
ture of about 2-3 MeV. As a consequence, it was in- 
ferred that pure hadronic stars are converted to strange 
stars or to hybrid stars within the first seconds after their 
birth. However, neutrino trapping in the protoneutron star 
phase strongly precludes the formation of a quark phase 
(Lugones and Benvenuto 1998, Benvenuto and Lugones 



1999, Vidaha, Bombaci & Parenti 2005). Then, it is pos- 
sible that the compact star survives the early stages of its 
evolution as a pure hadronic star. In this case, nucleation 
of quark matter would be triggered by quantum fluctua- 
tions in degenerate (T = 0) neutrino-free hadronic matter 
(see e.g. Grassi 1998, Iida and Sato 1998, Berezhiani et 
al. 2003, Harko, Cheng and Tang 2004, Bombaci, Parenti 
and Vidaha 2004, and references therein). 

A general feature of degenerate Fermi systems is that 
they become unstable if there exist any attractive interac- 
tion at the Fermi surface (Bardeen, Cooper and Schrieffer 
1957). This instability leads to the formation of a con- 
densate of Cooper pairs and the appearance of supercon- 
ductivity. In QCD any attractive quark-quark interaction 
will lead to pairing and color superconductivity (Barrois 
1977, Bailin and Love 1984 and references therein). Since 
the typical superconducting gaps in quark matter may be 
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as large as A ~ 100 MeV (see e.g. Afford 2001 and refer- 
ences therein) it is interesting to study the effects of color 
superconductivity in the process of nucleation. Some re- 
cent studies have gone in this direction (Drago, Lavagno & 
Pagliara 2004, Lugones and Bombaci 2005), but we shall 
give here a more complete and self-consistent study. 

In a recent work (Lugones & Bombaci 2005, here- 
after LB05) we have studied the formation of supercon- 
ducting quark matter in bulk, paying particular attention 
to the microscopic state of quark matter just after the 
deconfinement. As explained in LB05, several states are 
possible in principle. For example, quantum fluctuations 
could form a drop of /3-stable quark matter (hereafter the 
Q 13 phase). However, this is strongly suppressed with re- 
spect to the formation of a non /3-stable drop by a factor 
~ G For ' mi being N the number of particles in the criti- 
cal size quark drop. This is so because the formation of a 
/3-stable drop would involve the almost simultaneous con- 
version of ~ N/3 up and down quarks into strange quarks. 
For a critical size /3-stablc nugget at the center of a neu- 
tron star it is found N ~ 100—1000, and therefore the fac- 
tor is actually tiny. This is the same reason that impedes 
that an iron nucleus converts into a drop of strange quark 
matter, even in the case in which strange quark matter 
had a lower energy per baryon (Bodmer-Witten-Terazawa 
hypothesis). 

However, quantum fluctuations can form the so called 
Q*-phase bubbles, in which the flavor content of the quark 
phase is equal to that of the /3-stable hadronic system at 
the same pressure. Since no flavor conversion is involved, 
there are no suppressing Fermi factors, and a Q*-phase 
drop can be nucleated much easier. Once a critical size 
drop of the Q*-phase is formed the weak interactions will 
have enough time to act, changing the quark flavor frac- 
tion of the deconfined droplet to lower its energy, and a 
droplet of the (/3-stable) Q^-phase is formed. Notice that 
in degenerate matter, the intermediate phase can be made 
up by unpaired quark matter (hereafter the Q* unv phase) 
or by color-superconducting quark matter (Q* A phase). 
Which one of the two kind of droplets (Q1 np or Q*a) wm 
nucleate depends on the Gibbs free energy per baryon of 
each phase (g U np, 9 a), as discussed in detail in LB05. 

The analysis in LB05 has been made in bulk, i.e. with- 
out taking into account the energy cost due to finite size 
effects in creating a drop of deconfined quark matter in the 
hadronic environment. As a consequence of the surface ef- 
fects it is necessary to have an overpressure AP = P — Po 
with respect to the bulk transition point Po (see Fig. 1 in 
LB05). Thus, above Po, hadronic matter is in a metastable 
state, and the formation of a real drop of quark matter oc- 
curs via a quantum nucleation mechanism. In this work we 
shall study the quantum nucleation process including the 
effect of color superconductivity, and analyze the implica- 
tions for neutron stars. 1 

1 It is worth mentioning that Drago, Lavagno & Pagliara 
2004 have investigated the nucleation of a superconducting /te- 
stable phase made up by color-flavor-locked (CFL) quark mat- 



A final comment is worthwhile concerning the forma- 
tion of mixed hadron-quark phases, in which the electric 
charge is zero globally but not locally, i.e. the two phases 
have opposite charge (Glendenning 1992, 2001). As dis- 
cussed in Lugones and Bcnvenuto 1998, mixed phases can- 
not form in the here studied just- deconfined phase. The 
reason is that the flavor conservation condition guarantees 
that a just-deconfincd quark-matter drop has initially ex- 
actly the same electric charge than the hadronic drop from 
which it originated (i.e. zero). Of course, charge separa- 
tion could occur later on (if energetically preferred) and 
a mixed phase could form. However, notice that Debye 
screening effects and the surface tension can prevent mixed 
phases to form (see e.g. Tatsumi & Voskresensky 2003, 
Endo et al. 2005, Maruyama et al. 2006, and references 
therein). In any case, the study of mixed phases concerns 
the state of the system at times much larger than the ones 
that are addressed in this paper. 

2. The properties of the intermediate 
non-/3-stable quark phase 

As emphasized in the Introduction the intermediate phase 
can be in an unpaired state {Q* unp ) or in the so-called two- 
flavor color superconducting (2SC) state (hereafter Q* A ). 
In this section we resume the equation of state of the Q A - 
phase (for more details see LB05). 

The physical conditions imposed on the Q^-phase are 
flavor conservation and color neutrality of the quark gas. 
Flavor conservation means that the particle number (per 
baryon) of quarks u, d and s are the same in the hadronic 
phase and in the Q^-phase. The only difference is that 
in the hadronic phase quarks are confined inside hadrons, 
and in the Q* A -]A\a&e they are deconfined and paired. This 
condition can be expressed in terms of two parameters £ 
and n: 

n d =£,n u . (1) 

n s = n n u . (2) 

where is the particle number density of the i-species 
in the quark phase. The quantities £ = Yj 1 /Y^ and 
n = Y S H /Y^ are functions of the pressure, and charac- 
terize the composition of the hadronic phase. These ex- 
pressions are valid for any hadronic EOS. For hadronic 
matter containing n, p, A, E+, S°, £~, E _ , and 5°, we 
have 

^ _ n p + 2n n + n A + n s o + 2n s - + n E - 
2n p + n n + )ia + 2n s + + n s o + n s o ' 

n\ + n s + + n s o + n s - + 2n s o + 2n H - 

V = • (4) 

2n p + n n + n\ + 2n s+ + n s o + n s o 

Notice that £ and n determine univocally the number of 
electrons present in the system through electric charge 
neutrality of the deconfined phase: 

3n e = 2n u - n d - n s . (5) 

ter. However, as stated before, the direct nucleation of a drop 
of such a phase is disfavored by weak interactions. 
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We also impose two pairing conditions: n<j r 
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order to allow for paring between quarks d r with u g , and 
n ur = ndg in order to allow for paring between quarks 
u r with d g . Finally, the system must be globally colorless, 
that is: n r = n g = rib- 

From the above equations we obtain the number densi- 
ties of each quark species in the paired phase as functions 
of the flavor composition rj, £: 

4-2£ 
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I I , --ui ■ (9) 
The other particle densities are given by 

Tldg = TL uri and Tl S g = TL S r = ^sb- 

The pressure and Gibbs energy per baryon of the 
paired deconfined phase can also be written in terms of 
the same parameters: 
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where kj c = (f£t c — to 2 J 1 / 2 and \x = p ur . The chemical 
potentials (Xf c are obtained by inverting the following set 
of equations, 



n fc 



3vr 2 



2A 2 p 



fc ) 3 rfy , dg 



(13) 

(14) 
(15) 



fif c = (37r 2 n /c ) 1/3 / c = u b , d b 

(ifc = [(37r 2 n/ c ) 2/3 + m 2 s } 1/2 f c = s r , s 9 , s b 

where the number densities are given by Eqs. 

For sufficiently large values of A the energy cost in- 
vested in forcing the 2SC pairing pattern is compensated 
by the gain of the condensation energy and the pre- 
ferred phase is Q*^. For small values of A the preferred 
state is the Q* np -phase described in Olesen & Madsen 
1994, Lugones & Benvenuto 1998, and Bombaci, Parenti 
& Vidaha 2004. 

3. Quantum nucleation of quark matter bubbles 

To calculate the nucleation rate of quark matter in the 
hadronic medium we use the Lifshitz-Kagan quantum nu- 
cleation theory (Lifshitz & Kagan 1972) in the relativis- 
ts form given by Iida & Sato (1998). The QM droplet 
is supposed to be a sphere of radius 1Z and its quantum 
fluctuations are described by the lagrangian 



L(n,H) = -M{K)c 2 Jl - {K/c) 2 +M{1l)c 2 -U(K) ,(16) 



where A4(1Z) is the effective mass of the QM droplet, and 
U(1Z) its potential energy. 

Within the Lifshitz-Kagan quantum nucleation theory, 
one assumes that the phase boundary (i.e. the droplet sur- 
face) moves slowly compared to the high sound velocity of 
the medium (1Z << v s ~ c). Thus the number density of 
each phase adjust adiabatically to the fluctuations of the 
droplet radius, and the system retains pressure equilib- 
rium between the two phases. Thus, the droplet effective 
mass is given by (Lifshitz & Kagan 1972; Iida & Sato 1998) 



M(TZ) = 4ttph[1 - 



nb,Q* 
rihH 



11' 



(17) 



where pn is the hadronic mass density, rib t H and n^Q* are 
the baryonic number densities at a same pressure in the 
hadronic and Q*-phase, respectively. The difference in the 
Gibbs free energy is given by (Lifshitz & Kagan 1972, Iida 
& Sato 1998) 



U(TZ) = -nll 3 n K Q*(nQ* - fi H ) 



(18) 



where [in and /iq* are the hadronic and quark chemical 
potentials at a fixed pressure P. For comparison with pre- 
vious works notice that fj, is the same as the bulk Gibbs 
energy per baryon g = (p + P)/n B = niPi)/n B . For 
a = we have U(7l) = Volume xn B x Ag, and the bulk 
limit is recovered. 

Notice that we neglected the term associated with the 
curvature energy, and also terms connected with the elec- 
trostatic energy, since they are known to introduce small 
corrections (Iida & Sato 1998, Bombaci et. al 2004). The 
value of the surface tension a for the interface separating 
the quarks and hadrons phase is poorly known, and typical 
values used in the literature range within 10-50 MeV/fm 2 
(Heiselberg et al. 1993, Iida & Sato 1998). Notice that, 
since we are considering a as a parameter, color supercon- 
ductivity enters the above expressions only trough ptQ* 
and nb,Q* (see next section). 

The probability of formation of a bubble having a crit- 
ical radius can be computed using a semiclassical approx- 
imation. The procedure is rather straightforward. First 
one computes, using the well known Wentzel-Kramers- 
Brillouin (WKB) approximation, the ground state energy 
Eq and the oscillation frequency vq of the virtual QM drop 
in the potential well U(TZ). Then it is possible to calculate 
in a relativistic framework the probability of tunneling as 
(Iida & Sato 1998) 

A(E )- 



Po = exp 



where A is the action under the potential barrier 



(19) 



ME) 



{[2Mc 2 + E-U][U -E]} 1 ' 2 dTZ, (20) 



1Z± being the classical turning points. 
The nucleation time is then equal to 



(21) 
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Fig. 1. The critical mass of hadronic stars for the nu- 
cleation of a quark matter drop at the center. We have 
considered the parametrization GM1 (Glendenning & 
Moszkowski 1991) for the EOS of the hadronic phase, and 
two different values for the surface tension of the quark- 
hadron interface: a = 10 MeV/fm 2 and a = 30 MeV/fm 2 . 
Different curves correspond to different values of the Bag 
constant (indicated in units of MeVfm -3 ). For small val- 
ues of the pairing gap A the nucleated intermediate phase 
is unpaired (<2£„ p ) and is indicated with dotted line. For 
sufficiently large A the preferred intermediate phase is 
paired (Q%) and the critical mass decreases substantially 
(solid line). The maximum mass for hadronic stars with 
the GM1 EOS is MV}% X = 1.78M n . 



where N c is the number of virtual centers of droplet for- 
mation in the innermost region of the star. Following 
the simple estimate given in Iida & Sato (1998), we take 
N c = 10 48 . The uncertainty in the value of N c is expected 
to be within one or two orders of magnitude. In any case, 
all the qualitative features of our scenario will be not af- 
fected by the uncertainty in the value of N c . 

In order to explore the astrophysical implications of 
the nucleation we shall introduce the concept of critical 
mass of an hadronic star (see e.g. Bombaci et al. 2004). 
In bulk, the transition begins when the Gibbs conditions 
AP = and Afj, = are verified. However, as a conse- 
quence of the surface effects it is necessary now to have an 
overpressure AP = P — Pq > with respect to the bulk 
transition point. This overpressure will produce A/i < in 
Eq. {TBI allowing the barrier to be tunneled. Since there is 
an univocal relation between the mass of a compact star 
and its central pressure, we can consider that the nucle- 
ation time in Eq. (|21[) is a function of the mass of the 
star. The larger the overpressure, the easier will be to nu- 
cleate a bubble. In other words, the larger the mass of 
a pure hadronic star, the shorter the time to nucleate a 
quark drop at the center of the star. We define as critical 
mass M cr of the metastable hadronic star (HS), the value 
of the gravitational mass for which the nucleation time 



Fig. 2. The same as Fig. 1 but using GM3 for the hadronic 
phase. The maximum mass for hadronic stars with the 
GM3 EOS is MV}q x = 1.53M Q . 



is equal to one year: M cr = Mhs{ t — 1 Y r )- It is worth 
recalling that the nucleation time given in Eq. (|21|l is an 
extremely steeply function of the hadronic star mass (see 
e.g. Bombaci et al. 2004). Therefore, the exact value of the 
nucleation time (r = 1 yr in our case) chosen in the def- 
inition of M cr is not relevant 2 . For example, a variation 
of several orders of magnitude in the choice of r (e.g. 1 s 
or 10 17 s), will change M cr by less than 0.1 %. Therefore, 
pure hadronic stars with Mhs > M cr are very unlikely to 
be observed, while pure hadronic stars with Mhs < M cr 
are safe with respect to a sudden transition to quark mat- 
ter. M cr plays the role of an effective maximum mass for 
the hadronic branch of compact stars (see the discussion 
in Bombaci et. al 2004). While the Oppenheimer-Volkov 
maximum mass M^g x (Oppenheimer & Volkov 1939) is 
determined by the overall stiffness of the EOS for hadronic 
matter, the value of M cr will depend in addition on the 
properties of the intermediate non-/3 stable quark phase 
(Qunp or Q*a)- As emphasized in LB05, the critical mass 
exists even in the absence of surface effects, although its 
value is smaller. 

The results are shown in Figs. 1-4 where we show the 
critical mass of hadronic stars as a function of the differ- 
ent parameters of the equations of state. We have adopted 
rather common models for describing both the hadronic 
and the quark phase of dense matter. For the hadronic 
phase we used models which are based on a relativis- 
tic lagrangian of hadrons interacting via the exchange of 
sigma, rho and omega mesons. The parameters adopted 
are the standard ones (Glendenning & Moszkowski 1991). 
Hereafter we refer to this model as the GM equation of 

2 Notice also that, due to this tremendously steeply depen- 
dence of t with mass of the star, the poor knowledge of the 
factor N c in eq. (21) does not modify significantly the value 
obtained for M cr . 
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state (EOS) 3 . For quark matter we used the MIT Bag 
model with m u — m d = 0, m s = 150 MeV and a s — 4 . 

Before analyzing the results, it is worth to remark 
that there are two qualitatively different possibilities con- 
cerning the state of quark matter: either it is absolutely 
stable (i.e. it has an energy per baryon at P = and 
T = smaller than the mass of the neutron) or it is 
not absolutely stable. Since the actual case in nature is 
unknown the so called stability windows have been intro- 
duced (Farhi & Jaffe 1984). These windows show the re- 
gions in which each possibility is realized as a function 
of the different parameters of the quark model (e.g. the 
strange quark mass m s , the Bag constant B, and the pair- 
ing gap A). The stability windows for color- flavor locked 
quark matter have been presented in Lugones & Horvath 
(2002). Fixing the value of the strange quark mass to 
m s = 150 MeV it is found that paired /3-stable quark 
matter is absolutely stable for 



B < 73.1 MeV fm" 



l+,, ^ l ioow)' 



(22) 



The above expression is an approximation to order ml, 
accurate within a few percent in the range of interest (see 
Eq. (17) of Lugones & Horvath 2002). If quark matter is 
not absolutely stable, the stars containing quark matter 
are hybrid stars, i.e. containing /3-stable quark matter only 
at their interiors. If quark matter is absolutely stable, all 
the stars containing quark matter are strange stars, i.e. 
made up of /3-stable quark matter from the center to the 
surface. 

In Figs. 1 and 2 the effect of the pairing gap A is seen 
clearly. For small values of A the nucleated intermediate 
phase is in an unpaired state (dotted line). For sufficiently 
large values of A the gain of the condensation energy fa- 
vors the nucleation of a paired state (Q%) and the critical 
mass decreases substantially due to the lower energy cost 
of the phase. Notice that the effect of pairing is stronger 
for small values of the Bag constant B, where it can change 
M cr by more than a factor of 2. On the other hand, the 
variation of the critical mass due only to pairing effects 
(i.e. at constant B) is smaller than a ~ 20% in the regime 
of large B that corresponds to hybrid stars (i.e. non ab- 
solute stability of quark matter). 



3 The names GM1 and GM3 refer to the parameters given 
in the first and third lines of Table 2 of Glendenning & 
Moszkowski (1991). The corresponding compressibility and 
coupling constants are the following. For GM1: K — 300 MeV, 
(<? CT /m CT ) 2 = 11-79 fm 2 , (g„/m„f = 7.149 fm 2 , (g P /m p ) 2 = 



4.411 fm 2 , b = 0.002947, c 



-0.001070. For GM3: K 



240 MeV, {g„/m a y = 9.927 fm 2 , (g^/m^y = 4.820 fm 2 , 
(g P /m p f = 4.791 fm 2 , b = 0.008659, c = -0.002421. The 
two hadronic models include the lowest baryon octet (n, p, A, 
E+, E°, E-, E-, and S°). 

4 Other models, such as the Nambu - Jona - Lasinio model 
(see Buballa 2005), can give different results, specially in the 
low density regime for which the chiral symmetry is not com- 
pletely restored. 




a [MeV fm 



a [MeV frrT] 



Fig. 3. The critical mass of hadronic stars as a function of 
the surface tension a of the quark-hadron interphase. We 
show the results for two values of the Bag constant B (in- 
dicated in units of MeV fm -3 ). Dotted lines correspond to 
unpaired quark matter, solid lines to paired quark matter 
with A = 100 MeV. 



In figure 3 we see the effect of the surface tension a on 
the critical mass M cr . The effect of a is small for B larger 
than ~ 100 MeV fm -3 . Therefore, in this regime of B 
the bulk limit is a good approximation (c.f. LB05). For 
B smaller than ~ lOOMeV fm -3 the effects of <7 are im- 
portant, and increase their relevance as B decreases. This 
can be seen in Fig. 3 in the case for B — 60 MeV fm -3 , 
where M cr changes by a factor of 4 in the most plausi- 
ble range for a (between 10 and 30 MeV fm -2 ). However, 
notice that the curves tend to "saturate" for larger a. 

In Fig. 4 we explore the dependence of the results on 
the Bag constant B. We can identify three qualitatively 
different regions: 

(1) Region of Large B: For B larger than a critical 
value (~ 200 MeV fm -3 for GM1 and - 150 MeV fm -3 
for GM3) the critical mass M cr is no longer smaller than 
the maximum mass for hadronic stars M^g x . This im- 
plies that quark stars cannot form by means of quantum 
nucleation in this region of the parameter space. 

(2) Region of intermediate B corresponding to hybrid 
stars: For B between ~ 100 MeV fm -3 and ~ 150 - 
200 MeV fm -3 the critical mass M cr is very close (but 
smaller than) the maximum mass of hadronic stars. Notice 
that for this range of B, stars containing quark matter are 
hybrid since quark matter in /3-equilibrium has an energy 
per baryon (at P = 0) smaller than the neutron mass (we 
assume for simplicity that /3-stable quark matter is in a 
color flavor locked state). 

(3) Region of B corresponding to strange stars: For 
sufficiently small B quark matter in /3-equilibrium is ab- 
solutely stable (see Eq. (22)). As apparent from Fig. 4, the 
critical mass M cr is strongly dependent on B in this case. 
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Fig. 4. The critical mass of hadronic stars as a function 
of the Bag constant B. The results are shown for a = 
10 and 30 MeV fm -2 . Dotted lines correspond to un- 
paired quark matter, solid lines to paired quark matter 
with A = 100 MeV. The dashed horizontal line indicates 
the maximum mass of hadronic stars for each hadronic 
equation of state. 

4. Energy released in the conversion 

The amount of energy that can be released in the conver- 
sion of a pure hadronic star (HS) into a strange star or an 
hybrid star has been calculated in several previous works 
(Bombaci & Datta 2000, Berezhiani et al. 2003, Bombaci 
et al. 2004, Drago, Lavagno & Pagliara 2004). Moreover, 
Drago, Lavagno & Pagliara (2004) calculated the avail- 
able energy including the effect of color superconductiv- 
ity in the final star. However, since the critical masses in 
that work were not calculated employing the here stud- 
ied intermediate non-/3-stable phase, we present here new 
self-consistent results. 

As in previous works, we define AE as the difference 
between the gravitational mass of the initial hadronic star 
and that of the final (hybrid or strange) star having the 
same baryonic mass: 

AE = [M G (HS) - M G (QS)]c 2 . (23) 

According to the definition, AE is the total energy that 
can be released in the conversion. This energy can be lib- 
erated in an explosive manner (and give raise to GRBs) 
or in a less violent way (in which case observable signals 
would be less spectacular), but this depends on complex 
processes that are out of the scope of the present work 
(see e.g. Lugones, Ghezzi, de Gouveia Dal Pino & Horvath 
2002, Ouyed, Rapp and Vogt 2005, Paczynski and Hacnscl 
2005, Bhattacharyya, Ghosh and Raha 2005, and refer- 
ences therein). 

The results are shown in Fig. 5 and 6 where we show 
AE as a function of the paring gap A for different values 
of the Bag Constant B, using the equations of state GM1 



and GM3 for hadronic matter, and assuming that the /te- 
stable quark phase (for the final stellar configuration) is 
in the color flavor locked (CFL) state. We assume that a 
given hadronic star can convert into a quark star provided 
it has the critical mass. Those with a smaller mass cannot 
convert, and those with larger mass are already converted. 
Each point of the curves of Figs. 5 and 6 indicate the en- 
ergy released in the conversion of an hadronic star having 
the critical mass corresponding to that B and A. This 
hadronic star is converted into a quark star (having the 
same barionic mass) that can be a strange star or a hybrid 
star depending also on the values of B and A (this can be 
known with the help of Eq. 22). 

The released energy is mainly the result of two com- 
peting effects. For A above a threshold value (see Figs. 1 
and 2), pairing tends to diminish M cr , and therefore, de- 
creases the released energy since the conversion happens 
at the critical mass. On the other hand, there is an effect 
on the gravitational mass of the hybrid or strange star 
due to the condensation term in the equation of state of 
CFL quark matter. Depending on B and A the equation 
of state can be stiffer or softer than the unpaired case (see 
e.g. Lugones and Horvath 2003) and therefore AE can be 
either an increasing or a decreasing function of A. 

The dotted part of the curves in Figs. 5 and 6 corre- 
sponds to the case in which the intermediate non-/3-stable 
quark phase is unpaired. Therefore, the critical mass is in- 
dependent of A and the variations in AE are due only to 
the variation in the stiffness of the equations of state, that 
changes the gravitational mass of the final star. The full- 
line part of the curves corresponds to the case in which 
the intermediate phase is paired. Since the critical mass is 
a decreasing function of A (see Figs. 1 and 2) the released 
energy tends to decrease, specially for low values of B, for 
which the decrease of M cr is stronger (see e.g. the case 
of B = 60 MeVfm~ 3 ). For larger B the decrease in M cr 
is not enough to compensate the effect of the stiffness of 
the equations of state. In fact, for very large A, AE is al- 
ways an increasing function of A because the gain in the 
condensation energy compensates all other effects. 

5. Discussion and Conclusions 

In the present work we have studied the nucleation of 
quark matter drops at the center of cold deleptonized neu- 
tron stars. These drops can be made up by unpaired quark 
matter or by color superconducting quark matter. The na- 
ture of the nucleated phase depends on the details of the 
equations of state for quark and hadronic matter, and is 
relevant in the determination of the critical mass M cr of 
the hadronic stars above which a transition to quark mat- 
ter is possible. We have calculated M cr as a function of 
the several parameters of the quark model, and for differ- 
ent parametrization of the hadronic equations of state. In 
general, the dependence of M cr on B, A and a is mild if 
the parameters of the quark model correspond to hybrid 
stars, and strong if they correspond to strange stars. Also, 
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A [MeV] 

Fig. 5. The total energy released in the conversion of an 
hadronic star (having the critical mass) into a quark star 
(hybrid star or strange star depending on the values of 
B and A). In dotted lines we show the cases in which 
the intermediate non-/3-stable quark phase is unpaired, 
while the full lines corresponds to the case in which the 
intermediate phase is paired (see more details in the text). 
In this figure wc used the GM1 EOS for hadronic matter. 
The critical mass is calculated with a = 10 McVfm~ 2 . 
The Bag constant is indicated in units of MeVfm~ 3 . 




50 100 150 

A [MeV] 

Fig. 6. The same as Fig. 5 but employing the GM3 EOS 
for hadronic matter. 

the critical mass always decreases with A, and increases 
with B and a. 

As stated before, if quark matter is not absolutely sta- 
ble, the stars containing quark matter are hybrid stars, 
i.e. containing /3-stable quark matter only at their inte- 
riors. On the other hand, if quark matter is absolutely 
stable, all the stars containing quark matter are strange 
stars, i.e. made up of /3-stable quark matter from the cen- 
ter up to the surface. Moreover, it has been argued that if 
/3-stable quark matter is absolutely stable, then all neu- 



tron stars should be made up by quarks. The reason is 
that cosmic ray strangelets will be created after the merg- 
ing of a strange star with another compact star in a bi- 
nary system. The rate of binary mergers in the galaxy 
appears to be large enough so that even a conservative es- 
timate gives a large galactic production rate of strangelets 
( 10~ 10 Mo/yr), and correspondingly, a large flux in the 
interstellar medium. When a strangelet penetrates a neu- 
tron star, it will grow by absorbing free neutrons, convert- 
ing the whole star into strange quark matter. Similarly, all 
massive stars that accumulate a strangelet in its core will 
give birth to a strange star rather than a neutron star af- 
ter the explosion as a core collapse supernova. Because of 
these arguments, it has been argued that if some neutron 
stars are undoubtedly identified as being conventional, 
the inevitable conclusion is that strange quark matter is 
not absolutely stable and there exist only conventional 
neutron stars or (eventually) hybrid stars (see Caldwell 
& Friedman 1991, Balberg 2005). However, notice that 
it may happen that the strangelet contamination of the 
galaxy is suppressed, for example, if strange quark stars 
are disfavored in binaries (Bclczynski et al. 2002 and ref- 
erences therein). On the other hand, strangelets could dis- 
integrate when impacting onto the external layers of the 
star. But, the physics of strangelet fragmentation is largely 
unknown and it is not clear whether it is effective for all 
incident energies (specially for those strangelets impacting 
with low energy). In conclusion, if quark matter is in fact 
absolutely stable, it is not clear whether the conversion 
would be triggered from inside (quantum nucleation) or 
from outside (strangelet contamination), and therefore, it 
is uncertain whether the critical mass studied in this paper 
is meaningful in the corresponding region of the parameter 
space (B, A, etc.). 

On the other hand, if quark matter is not absolutely 
stable, strangelets do not exist. Therefore, the trigger of 
the conversion to quark matter cannot come from out- 
side the neutron star, and the here studied mechanism is 
the most plausible one to trigger the conversion in cold 
deleleptonizcd neutron stars. In this context, a relevant 
result presented here is that for a large region of the pa- 
rameter space corresponding to hybrid stars, the critical 
mass is very near (but smaller than) the maximum mass 
of hadronic stars. This means that all compact stars with 
masses up to M cr should be normal hadronic stars while 
those with masses above M cr may be hybrid stars (if they 
do not collapse into a black hole). In the rest of the pa- 
rameter space corresponding to hybrid stars the critical 
mass M cr is no longer smaller than the maximum mass 
for hadronic stars M^g x . This implies that quark stars 
cannot form by means of quantum nucleation in this re- 
gion of the parameter space, and probably they don't exist 
at all (if other mechanism of quark matter formation does 
not operate). 

We have also calculated the amount of energy that 
can be released in the conversion of an hadronic star into 
a quark or hybrid star (assuming that the /3-stable quark 
matter phase is in a color flavor locked state). The total 
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released energy is in the range 3 x 10 52 erg - 4 x 10 53 erg, 
and therefore it is sufficient to power a gamma ray burst. 
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